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Report  Title 

Organizing  Principle  for  Tipping  Points  in  Social  Networks 

ABSTRACT 

IN  this  talk  ppt  given  at  Dr  Mubarak  Shah’s  Center  for  Research  in  Vision  at  U  Central  Florida,  Orlando,  March  25  - 
26  2013,  we  describe  the  overarching  principles  that  guides  the  recent  work  on  Tipping  points  and  committed 
Minorities  in  Signalling  Multi-Agent  Social  Networks 

We  find  that  a  scalar  stochastic  Differential  equation  can  be  derived  to  give  good  estimates  of  the  first  exit  times  such 
as  consensus  times  and  to  study  the  influence  of  diehards. 

For  very  large  populations  we  find  that  a  one-dimensional  center  manifold  on  which  the  mean-field  dynamics  in  slow 
time  is  easily  shown  to  consist  of  nodes  and  saddle  and  heteroclinic  orbits  linking  these  equilibria  -  the  saddle  node 
bifurcation  then  led  to  an  easy  determination  of  the  critical  value  of  diehard  fraction  needed  to  accelerate  the  network 
to  consensus  in  the  minority  opinion. 

For  smaller  crowds  on  the  order  of  1000  -  5000  agents,  we  find  that  the  demographic  noise  reflected  in  simulations 
by  the  significant  variance  in  consensus  times,  must  be  taken  into  account  to  determine  accurate  estimates  of  the 
expected  value  of  consensus  times  and  to  study  the  effects  of  committed  minority. 
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Tipping  Point  of  NG 

A  minority  of  committed  agents  can  persuade 
the  whole  network  to  a  global  consensus. 


The  critical  value  for  phase  transition  is  called 


Saddle  node  bifurcation 


Below  Critical 


Above  Critical 


Meanfield  Assumption  and  Complete 
Network 


The  network  structure  is  ignored.  Every 
node  is  only  affected  by  the  meanfield. 

The  meanfield  depends  only  on  the 
fractions(or  numbers)  of  all  types  of  nodes. 

Describe  the  dynamics  by  an  equation  of 
the  meanfield  (macrostate). 


Scale  of  consensus  time  on  complet 
graph 


2  word  Naming  Game  on  complete  graph 


Expected  Time  Spend  on  Each 
Macrostate  before  Consensus  (without 

committed  agents) 


n. 


n 


B 


0  0 


NG  with  Committed  Agents 

'1  q=0.06<q<,  q=0.12>qc 


q  is  the  fraction  of  agents  committed  in  A. 

When  q  is  below  a  critical  value  q,,,  the  process 
may  stuck  in  a  meta-stable  state  for  a  very 
long  time. 


2  Word  Naming  Game  as  a 
2D  random  walk 


P{A-)  =  UAiN  -nA+  nB)/2N^ 
n_B  P{A+)  =  {N -UA- nB){N +  nA- nB)/2N‘^ 

P{B-)  =nB{N  +  nA-  nB)/2N'^ 

P{B+)  =  {N-nA-  nB){N  -uaA  nB)/2N^ 

Po  —  {riA  +  nB)/2N  +  {ua  -  ubY /2N‘^ 


P(B+) 


0 

Transient  State 

-> 

• 

Absorbing  State 

n  A 

Linear  Solver  for  2-Name  NG 


absorbing  time  starting  from  network  state  {nA-Jin)  ■ 
t{nA-Jf'B)  is  the  average  time  stay  in  network  state  before  leaving. 

Have  equations: 

P  (.4+  )t  (ri^  +  1 ,  fiB )  +  ^ (-4- )r  -  1 ,  )  +  P (B + )r  (r^ ,  +  1 )  +  P (5- )t (tia,  tib  -  1 ) 

r[nA,riB}  =  - j— - 

Then  we  assign  an  order  to  the  coordinates,  make  r{riA,nB)  , 
into  vectors,  and  finally  write  equations  in  the  linear  system  form: 


r  =  Mr  +  t 


SDE  models  for  NG, 


dX^  =  fddt  -h 


NG  and  NG 


Diffusion  vs  Drift 


Diffusion  scales  are  clear  from  broadening 
of  trajectories  bundles 

Drift  governed  by  mean  field  nonlinear 
DDEs  can  be  seen  from  the  average  / 
midlines  of  bundles 


Assume  a  very  natural  social  -  political  condition,  generalizing  NG 

Where  the  network  is  divided  into  k+1  sub-populations 

Each  with  a  different  fixed  propensity  to  signal/vote/utter  the  opinion  A 

When  hearing  the  word  A  a  node  from  subgroup  j  <  k  will  move  to  subgroup  j+1 
Likewise  hearing  B  a  node  from  j  >  0  will  move  to  subgroup  j-1 

The  probability  that  a  node  from  subgroup  j  will  signal  A  is  j/k 

Same  node  has  probability  1  -  j/k  of  signalling  word  B 


Additional  rules  of  k-NG 


A  node  s  is  drawn  at  random  from  the 
network  and  sends  out  a  signal  A  with  p 

or  s  signals  B  with  prob  =  1  -  p 

Next  a  node  L  is  drawn  at  random  to 
receive  the  signal  A 

or  node  L  is  drawn  to  receive  the  signal  B 


k-NG 


The  subgroups]  =  0,  k  correspond  to 
those  nodes  that  are  completely  convinced 
of  the  B  and  A  opinion  resp.  OR 

Equivalently  those  nodes  that  with  prob  1 
Signals  B,  A  resp. 


Voting  or  Polling 


The  key  network  quantity  is  an  average 
network  opinion  obtained  by  polling: 

p  =  sum  over  subgrps  j  =  0  to  k 
of  n(j)j/kN 

It  gives  the  probability  of  a  speaker  chosen 
at  random  signalling  the  word  A 


Stochastic  Dynamics 

Derivation  of  k-dim  coupled  Random  Walk: 


For  j  =  1  to  k-1, 

nC.t+1)  =  n(j,t)  +  1,-1  with  resp.  prob. 
P(+1)  =  p(t)  nC-1,t)/N  +  (1-p(t))  n(j+1,t)/N 
P(-1)=  (1-p(t))nC,t)/N  +  p(t)nC,t)/N 
=  nC,t)/N 


Random  walk 


For  the  distinguished  subgrps  at  both  ends 
of  the  opinion  spectrum,  j  =  0  ,  k 

n(0,t+1)  =  n(0,t)  +  1, 0,  -1  with  prob. 
P(+1)  =  (1-p(t))  n(1,t)/N 
P(0)  =  (1-p(t))  n(0,t)/N 
P(-1)=  p(t)n(0,t)/N 


Subgroup  k 

n(k,t+1 )  =  n(k,t)  +1,0, 

P(+1)  =  p(t)  n(k-1,t)/N 
P(0)  =  p(t)  n(k,t)/N 
P(-1)  =  (1-p(t))n(k,t)/N 


-1  with  prob. 


Shadow  Walk 


Define  a  scalar  nonlazy  random  walk  on 
the  quantity  p  whose  consensus  times  are 
lower  bounds  for  the  actual  expected 
times. 

p(t+1)  =  p(t)  +  1/kN,  -1/kN  with  prob 


P(+1)  =  p(t) 
P(-1)=  1-p(t) 


SDE  or  Diffusion  Model 


Ill  ooiitiijiioiis  time  rlu’  SDE  is 


M  =  \2P,-[)/k 


Solution  of  SDE 


now  using  stjuidanl  Kolmogorov  BiMtkwiuds  ['([Uiirioii  iiiHliod 


V  Ohi  Oti 

Tv 

(}.«((,?)  =  =  Pr{/'  = 


iiip)  =  linn((;>,/) 

f—  It 


Expected  Times  to  Consensus 


Fii’st  calculate 


=  exp  [  2  v'iV 


Next  we  get 


/(f  _  io  exp  [2^/^Wx  -  dx 

foG{x)dx  exp  (2\f ^W ^  -  ^^^dx 


Expected  times  -  exit  times  -  stop  times 


and  the  expected  time  to  consensus  of  the  A  o]>iiiion  without  committed 
agents  and  conditioned  on  eventual  fixarioii  of  rlic  A  opinion  is  given  hy 


\\\mv 


is  tlie  unconditioned  expected  time  to  consensus  given  tliat  Po  =  p  and  solves 
tlie  stationaiv  form  of  the  KBE 


(PT  2M  dT  2u(p) 

^\ith  boimdan'  conditions 

•h 


liiM 

p— n 


t  <  X 

0. 


continued 


The  solution  is  given  by 


/ 


\ 


rlj  CG(^)'i^l'G(x)dr 

gG(xMxJ,rfv'" 


Higher  stubbornness  -  same  qualitative, 
robust  result 
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Higher  stubbornness  -  same  qualitative, 
robust  result 
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other  NG  variants  -  same  1 D  manifold 


Nq 


Figure  2:  Si]iiula,tioii:^  of  tiie  Naming  Game  with  K=1U  and  N  =  5UIJ.  Some  tra- 
jtv'TfjrH'H  (ici  not  {4 [>i:iroa4'li  tli<'  rfiiH'i'  m^Liiafolc]  ;ih  cj^iiiclcly  aw  l  (3i<l  in  tb(’ 
case.  Those  simulations  that  start  with  p=0.5  approach  the  centor  manifold  he- 
loi'e  tlriftiiig  to  the  consensus  state.  Howevner.  the  aiinulat ioiLs  tliat  start  w^aili  a 
less  cent i  nl iznd  p  value  rlrift  to  tlie  consensus  state  more  quickly.  Note  tlial  tiic 
simulations  starting  in  the  polai'ized  s^rate  with  p=0.4  come  clo^  to  the  center 
luaiill^ j|<.]  1jch>rc'  (Irirtiii^  to  cx iliseiLsiLs.  On  1  lif?  ot  her  IijiikI.  llii^  sjiiiM LL[i.t iems  iu 
cent  rallied  estate  with  p=0.'6  drift  to  coiiseixsns  hefore  any  significant  number  of 
agentsi  enter  A'q. 


3D  plot  of  trajectory  bundles  - 
stubbornness  K  =  10  as  example  of 

_ _ ? _ X  /XX  ^ _ ?x _ _  _ _i  I  ? _  r\r\Jk  r\\ 


K=  ia.N  =  5(» 


urtfom,  Pj  =  D  5 
centered.  Pj  = 
poanzed,  =  0.5 
centered,  Pj  =  0.6 
po-anzed,  p^  =  0.4 


Fi^iirr  4:  A  3D  plot  of  tho  siiiuilanoiis  nf  the  Naming  Ganip  T.\-nh  K=3  and  N 
=  Hf'rc,  all  Hiimibuii approach  rlie  cnHcr  maiiifold  licfore  p  de\i:Uf'H  too 

much  hoiti  its  initial  vnluc. 


Consensus  time  distribution 


Recursive  relationship  of  P(X,  T),  the  probability 
for  consensus  at  T  starting  from  X,  Q  is  the 
transition  matrix. 


P{nAynB,T  +  1)  =  Q(?i^  +  l,nB\nAynB)P{nA  +  l,nB,T) Q(71a  - 

-h  l^T)  +  -  1,1) 

+Q(7M,?^b  +  2|7i4,nB)P(n4,7iB  -\-2,T)  +  Q{nA  +  2,77BkM,riB)P'(ri^  +  2,71^,!) 

Take  each  column  for  the  same  T  as  a  vector: 

P{T  +  1)  =  Q*P{T) 

Calculate  the  whole  table  P(X,T)  iteratively. 

Take  each  row  for  the  same  X  as  a  vector: 

^{nA,nB)i'^c  =  T)  =  {P{nA,nB,T)). 


Consensus  time  distribution 


"^c  X  10“*  X  10’’ 


Red  lines  are  calculated  through  the  recursive  equation. 

Blue  lines  are  statistics  of  consensus  times  from  numerical  simulation(very  expensive), 
(done  by  Jerry  Xie) 


P(T)  *  std(T) 


p=0.12 


Consensus  Time  distribution 


p=0.08 


(  T -E[TJ  )  /  std(T^) 

Below  critical,  consensus  time  distribution  tends  to  exponential. 
Above  critical,  consensus  time  distribution  tends  to  Gaussian. 

For  large  enough  system,  only  the  mean  and  the  variance  of  the 
consensus  time  is  needed. 


NG  on  RGG 
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travelling  wave  solution  of  u(x)  when  q=0.05<q^ 
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NG  on  RGG  past  Tipping  point 
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Homogeneous  Pairwis  Assumption 

The  mean  field  is  not  uniform  but  varies  for  the  nodes 
with  different  opinion. 


Make  it  rigorous: 

suppose  three  nodes  are  linked  as  1-2-3 
Xi  is  the  opinion  of  node  i 
ki  is  the  degree  of  node  i 
P{X,  X2)  ^  P{Xi) 

P{X,X2,X^)  =  P{X,X2] 

£[ftllxl]=<  k  > 

P{X,\k,)  =  P{X,) 

P{X,X,,k^,k2)  =  P{X,X2) 


Numerical  comparison 


Trajectories  mapped  to  2D  macrostate 


0.95 


Concentration  of  the  consensus  time 


Change  of  the  tipping  point  w.r.t.  the 
average  degree 


(hb 


N=10000 


*  MeanField 
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The  local  mean  field  for  the  node  with  opinion  C\ 

P(  |C)  =  [P{A\C),P{B\C),P{AB\C)f  ,C  ^  A,B,AB 
The  number  of  different  type  of  links: 


(  P(^|A)  \ 
^(L)  =  F(B|A) 
VT(A5|A)/ 


1 

2La-a  La -b  +  La-ab 


2ivl-v4\ 
La-b  ] 
La-ab ) 


P{\B){L) 


f  P(AB)\ 
P(BB) 
\P(ABB)J 


1 

La-b  +  +  Lb-ab 


La-b  \ 

'2Lb-b  1 


F(.|AB)(L) 


/  P(A|A5)  \ 
PiBAB) 
\P{AB\AB)J 


1 

La-ab a  Lb- ab  a 2Lab-ab 


(La-ab 
Lb-ab 
2Lab-ab 


Analyze  the  dynamics 


1  .Choosing  one  type  of  links,  say  A-B,  and  A  is  the  listener. 

2.  Direct  change:  A-B  changes  into  AB-B. 

3.  Related  changes:  since  A  changes  into  AB,  <k>-1  related 
links  C-A  change  into  C-AB.  The  probability  distribution  of  C 
is  the  local  mean  field  P(  |A). 


Direct 


Local  mean  field  equation 


£[AL|L]  =  ^[D  +  (<  k  >  -1)R]  L 
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R  =  ( 0,^[g^P(.|A)  +  QbP(-\BIQa\^P{^)  -  -/{-  ABUQBljPm  -  -^P(-\ABI-{Qa  +  Qb)P(-\AB)  ) 
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Normalized  equation: 
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<k> 


<  k  >  —1 
<k> 


<k> 
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SDE  model  of  NG 


r+AT-1 

AX{AT)  =  ^ 


i=T 


£J[A?i] 


S(An) 


dX  =  fidt  +  -^dW 


[dX,  dX] 


N 


Merits  of  SDE  model 


dX  =  fidt  +  -^dW 

v"iV 

Include  all  types  of  NG  and  other 
communication  models  in  one  framework 
and  distinguish  them  by  two  parameters. 

Present  the  effect  of  system  size  explicitly. 

Collapse  complicated  dynamics  into  1-d 
SDE  equation  on  the  center  manifold. 


Thanks 


